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ABSTRACT 

O 
O 

CN ' We present high resolution interferometric and single dish observations of 

molecular gas in the Serpens cluster-forming core. Star formation does not 

appear to be homogeneous throughout the core, but is localised in spatially- 

and kinematically-separated sub-clusters. The stellar (or proto-stellar) density 

■ in each of the sub-clusters is much higher than the mean for the entire Serpens 

cluster. This is the first observational evidence for the hierarchical fragmentation 

> ■ of proto-cluster cores suggested by cluster formation models. 
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1. Introduction 



It is generally accepted that most stars are born in clusters (cf. Clarke et al. |2000| ). 



The way in which clusters form and evolve is therefore likely to influence the distribution 



of masses for stars in the field, the initial mass function, IMF (Salpeter |1955| ; Scalo |1986| ) 



In very young, embedded clusters the distribution of stellar masses is often similar to the 
IMF (Palla & Stahler 1999j ; Meyer et al. |2(J(J(J| ), and the mass spectra of prestellar and 



protostellar condensations in the Serpens and p-Ophiuchi cluster-forming cores are also 
consistent with the IMF (Testi & Sargent [199^ , hereafter TS98; Motte et al. [19981) . 

It has been suggested that within stellar clusters star formation occurs preferentially 
in sub-clusters where the stellar density is much enhanced (Clarke et al. |2000| ). This 
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has important implications for cluster evolution. For example, more massive stars could 
be produced by coalescence (Stahler et al. |2000| ). To date, there is little evidence for 
sub-clustering in the Orion Nebula Cluster (Bate et al. |1998j ) or in the smaller clusters 
around intermediate mass pre- main sequence stars (Testi et al. |1999 ). However, models 
advocate sub-clusters with a much higher stellar density at the time of formation (Bonnell 



et al. |1998| ). It is therefore important to establish if sub- clustering is present in the very 
youngest clusters. The mean stellar densities and the stellar to gas mass ratio in such 
sub-clusters can provide critical observational constraints on coalescence, competitive 
accretion and binary evolution models (Bonnell et al. [1997]; |1998| ; Kroupa |1999| ). 



The Serpens molecular cloud, at ~ 310 pc (de Lara et al. |1991| ), is one of the most 
active nearby cluster-forming cores. Inside the 500-1500 M , ~0.6 pc diameter cloud of 
molecular gas (White et al. |1995| ) is a young protocluster comprising about one hundred 
embedded young stellar objects (YSOs), protostars and prestellar clumps (Strom et 
al. [T976| ; Eiroa & Casali [T992| ; Giovannetti et al. p98| ; Kaas |1999| ; Casali et al. [1991 ; Hurt 
Sz Barsony |1993| ; TS98). Numerous jets and molecular outflows have also been detected 
(Rodriguez et al. |1989| ; White et al [1991 ; Eiroa et al. [T997| ; Herbst et al. |1997| ; Wolf-Chase 
et al. |1998| ; Davis et al. |1999| ; Hodapp |1999| ; Hogerheijde et al. |1999| ). The total estimated 
mass of the YSOs, protostars, and prestellar clumps is in the range 40-80 M Q (Giovannetti 



et al. 1998 ; TS98), implying an overall star formation efficiency of 2-5%, similar to most 
nearby molecular clouds. The proto-cluster radius, ~ 0.2 pc, and mean stellar density, 
~ 400-800 stars/pc 3 , are typical of very young embedded clusters (Testi et al. |1999| ), making 
this an ideal laboratory for studying early cluster formation processes. Here, we present 
wide field, high resolution, aperture synthesis and single dish millimeter-wave molecular 
line observations of the Serpens core which support the concept of at-birth sub-clustering. 



2. Observations and results 

Owens Valley Radio Observatory (OVRO) millimeter wave array mosaic observations 
of the 5.5' x 5.5' inner region of the Serpens molecular core in the CS(2-1) transition at 
97.98 GHz were obtained at the same time as the 3 mm continuum data described by 
TS98. Details of the observations, data reduction, and mosaic construction are described 
elsewhere (TS98; Testi & Sargent |2000| ) . Spectral resolution was 0.4 km/s over a 24 km/s 
bandwidth centered at Vlsr =8 km/s. The synthesized beam is 5".5 x 4".3 (FWHM), 
and the noise level in each channel of the final cleaned mosaic cube is ~ 140 mJy/beam. 
During spring/fall 1998, additional observations of SMM4, S68N and the CS6, CS7 and CS8 
features identified in Figure [lja were acquired (Testi & Sargent |2000|) . The digital correlator 
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was configured so that the N 2 H + (l-0) transition at 93.2 GHz and the CH 3 OH(2_ 1 -l_ 1 )E2 
and (2 -lo)A + transitions at 96.7 GHz were detected simultaneously in the lower and upper 
sidebands respectively. Spectral resolution was 0.4 km/s over a 25 km/s band. Calibration 
was carried out as in TS98. The datasets were imaged using the AIPS IMAGR task. 

A 10' x 10' region of the Serpens cloud core, centered on a = 18 h 27 m 20 s , 5 = ri2'30" 
(B1950.0), was also mapped in the N 2 H + (l-0) line in March and October 1999, using 
the 13.7-m telescope of the Five College Radio Astronomy Observatory^ (FCRAO) and 
the SEQUOIA focal plane array. For the one beam sampled maps, the half power beam 
width (HPBW) was ~52". Low noise InP MMIC-based amplifiers resulted in mean receiver 
temperature of 70 K (SSB), and system temperatures of 160-250 K. The spectrometer was 
an autocorrelator with 24 kHz (0.077 km/s) spectral resolution and 20 MHz bandwidth. 
Typical integration times were 5 to 15 minutes in frequency switching mode, with a throw 
of 8 MHz. The main beam efficiency, to convert the antenna temperature to brightness 
temperature, ?? m b, is 0.51, and the final rms was ~0.05 K (T m &). 



2.1. OVRO maps 

In Figures |l]a and |I|b we show the CS(2-1) integrated intensity mosaic overlaid on our 
3mm continuum map (TS98) and the CS(2-1) mean velocity mosaic, respectively. In Fig. |l|a 
some of the brightest CS(2-1) features are labelled to simplify discussion; the location of 
the NIR cluster surrounding SVS-2 and SVS-20 (Strom et al. |1976j ) is also indicated. In 
Fig. |l]b, the orientations of the jets emanating from SMM1, SMM3, SMM4, and A3, as well 
as the NIR reflection nebula associated with SMM5 are marked. 

At most, 20% of the flux detected in single dish maps of the optically thick and 
extended CS(2-1) emission (cf. McMullin et al. |1994j ) is recovered. There is therefore little 



point in comparing the CS and the optically thin 3 mm emission (TS98). However, the 
molecular line observations provide information about the outflows emanating from the 
various millimeter sources; the abundance of CS can increase by a factor of almost 100 
in outflows (Bachiller & Perez Gutierrez |1997| ), and most of the clumps in Figure [l]a are 



likely to be compact enhancements in shocked regions. In most cases the CS linewidths 
are relatively broad, >2-3 km/s, and mean velocities differ by 3-4 km/s from the systemic 
velocity of the Serpens core (Figure [I]b). Strong and broad CH 3 OH emission, typical of 
shocked material, was also seen in our followup observations (Figure 2), and coincides with 
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the CS features. By contrast N 2 H + emission is spatially coincident with the continuum 
sources, as expected for a tracer of cores and envelopes. 



2.2. FCRAO N 2 H+ map 

Figure 3 shows the FCRAO/SEQUOIA large scale channel maps of the N 2 H + (l i-0i 2 ) 
isolated hyperfine component. Four cores, labelled A, B, C, and D in order of increasing 
Vlsr, and an extended, lower surface brightness "spur" can be identified. The spur 
stretches from the NW (C, D) to SE (A, B) cores and continues beyond these. Cores A and 
B are spatially associated with the SMM3/SMM4 region, while C and D encompass the 
SMM1/S68N region. Spectra at the emission peaks in A, B, C, and D are also presented in 
Fig. 3. These positions are given in Table 1 as offsets from a = 18 h 27 m 20 s , 5 = 1°12'30". 
Also listed are the Vlsr, hne width, AV, and virial mass, M^, for each core. The fraction 
of the total N 2 H + flux observed with the array is iovRO', is the total mass of 3 mm 
continuum cores within each gaseous core (TS98), and M* is the total mass of YSOs, 
derived assuming a mean stellar mass of 0.3 M (Giovannetti et al. |1998| ; Kaas |1999| ). The 



FCRAO and OVRO values of AV and Vlsr are consistent, suggesting that the array is 
detecting a fraction of the emission from extended envelopes, rather than the compact cores 
found by TS98. 



2.3. Identification of outflow sources 

The powering sources of the CS outflows were identified from published optical, 
infrared and millimeter observations, and are represented by open diamonds in Figure 1. 
Features CS3 and CS4 as well as CH 3 OH emission (Fig. 2) suggest a compact outflow at 



p.a.~140°, centered on the S68N continuum source (cf. Wolf-Chase et al. |1998|) . CS2 could 
either be an extension of the S68N flow or the counterflow of the H 2 jet from mm source 
A3 (Hodapp |1999| ; TS98). Based on its alignment with an H 2 jet (Hodapp |1999|) and the 



locations of HH460 and a CO(2-l) bullet (Davis et al. |1999| ), we associate CS1 with an 
outflow from SMM1. CS6, CS7, and CS8 probably mark the interaction of the counterflow 
from SMM1 with the B molecular core. All these outflows are oriented along p.a.~140°. A 
double-peaked CS component reported by Williams & Myers ( |1999| ) at the position of CS1 
was not detected; it is probably extended and resolved out in our observations. 

In the south-east, CS9 and CS10 are oriented like CH 3 OH (Fig. 2) and H 2 knots near 
SMM4 (Eiroa et al. |1997| ) along p.a.~180°, and are probably part of an outflow centered 
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on that source (c.f. Hogerheijde et al. |1999|) . Likewise CS12 is elongated along p.a.~170° 
chain of H 2 knots associated with SMM3 (Herbst et al. |1997| ). Davis et al. ( |1999[ ) 
noted that two Herbig-Haro objects and SMM3 are aligned along a position angle which is 
almost orthogonal to the chain of H2 knots. However, the kinematic properties of the knots 
suggest no physical connection with the HH objects. 



3. Discussion 

Our new interferometer and single-dish maps of the Serpens core indicate sub-clustering 
at an early epoch of cluster formation. Three separate properties argue for sub-clustering, 
spatial segregation, outflow orientations, and circumcluster gas kinematics. Approximately 
one third of the near infrared cluster members are concentrated in a ~0.1 pc radius region 
surrounding SVS-2 and SVS-20 (Giovannetti et al. |1998| ; Kaas |1999| ), while the millimeter 



and sub-millimeter sources are largely concentrated in the SE (A, B) and NW (C, D) 
fragments (TS98; Davis et al. |1999|) . The orientations of the outflows observed in the SE 
and NW are quite different. In the NW, all three flows are oriented along p.a.~140°, as 
is the near infrared reflection nebula centered on SMM5 (Kaas |1999| ; see also Figure 0b). 



In the SE, the two outflows are aligned approximately north-south, with mean p.a.~175°. 
In Figure 3, the NW and SE sub-clusters are embedded in discrete N 2 H + (l-0) clumps, 
separated in velocity by ~1 km/s (Table 1). Each clump comprises two cores: peak 
velocities of cores A and B differ by only 0.5 km/s (~ 1/2 AV), C and D by 0.1 km/s 
(~ 1/10 AV). Thus the spatially distinct clumps are also kinematically separated, while 
there is reasonable internal velocity coherence. 

It appears that the Serpens core encompasses at least three sub-structures - the 
NW and SE sub-clusters and the NIR cluster. Star formation is currently occurring 
simultaneously in the NW and SE sub-clusters, both of which contain roughly equal 
fractions of prestellar, protostellar and infrared sources. The NIR cluster is probably more 
evolved but we see no evidence of the progressive pattern of star formation proposed by 
Casali et al. ( |1993j ). The kinematics and outflow orientations indicate that each subcluster 



originated in a separate fragment of the cloud which subsequently fragmented into the 
smaller cores seen in the N 2 H + (l-0) maps. Within these are the 3 mm continuum cores that 
are likely progenitors of single stellar systems (TS98). Taken together these observations are 
consistent with either the hierarchical fragmentation picture advocated by Elmegreen ( |1997 ; 



1999| ) or the spontaneous fragmentation suggested by Myers ( |1998j ) for the formation of 
stellar clusters, and provide the first observational support for hierarchical fragmentation 
within a cluster-forming core. We note that while the mean stellar /protostellar density of 
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the entire Serpens core is ~400-800 stars/pc 3 , most of the proto-cluster members are within 
the three subclusters where densities reach 2000-4000 stars/pc 3 . It is very likely that, 
within a few million years, the cluster will evolve to a size and mean density very similar to 
those of embedded clusters around intermediate mass stars. 
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Table 1. N 2 H+(l-0) Cores Parameters 
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Fig. 1. — a) Integrated intensity contours for CS(2-1) emission in the Serpens core are 
shown in green, overlaid on the 3 mm continuum map (TS98). Contour levels begin 
at 0.9 Jykm/s/beam (3<r) and are spaced by 0.3 Jykm/s/beam to 3 Jykm/s/beam and 
thereafter by 1 Jykm/s/beam. The dashed ellipse represents the NIR cluster (Kaas |1999| ). 
The CS features discussed in the text are marked CS1 to CS12. b) Color-coded map of 
the mean velocity variations across the CS structures of Fig. |l|a. Known outflow sources, 
represented by open diamonds, are labelled and solid lines indicate the orientations of 
observed jets or, for SMM5, reflection nebulosity. The dashed line from SMM1 indicates 
the direction to HH460. The two rectangles mark the areas shown in Figure 2. Increased 
noise in the upper left and lower right corners of both images is due to loss of sensitivity at 
the edge of the mosaic. 
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Fig. 2.- Integrated intensity maps of CH 3 OH(2 fc -l fc ) and N 2 H + (l-0) for S68N and 
SMM4 the two regions indicated in Figure |l]b. Synthesised beams (FWHM) are shown 
by black ellipses in the lower right corners. Contour levels start at 3a and are spaced by lcr 
(0.3 Jykm/s/beam for S68N and 0.2 Jykm/s/beam for SMM4). 
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Fig. 3.- Top panel: FCRAO/SEQUOIA N 2 H + (l 01 -0i 2 ) channel maps. The five main 
kinematical structures, A, B, C, D, and the "Spur" are each labelled at their approximate 
Vlsr, which is given in the top-left corner of each panel. Contour levels are: 0.18 to 2.18 by 
0.36 K. The square in the bottom left channel shows the area mapped at OVRO. Bottom 
panel: N 2 H + (l 01 -0i 2 ) spectra observed towards the four cores A, B, C, and D, as defined 
in Fig. 3. In each panel the offset in arcsec from position a = 18 h 27 m 20 s , 5 = 1°12'30" 
(B1950.0) is given. 



